Transient response of a quantum wave to an instantaneous potential step switching 



F. Delgado, 1 H. Cruz, 2 and J. G. Muga 1 

1 Departamento de Quimica-Fisica, Universidad del Pais Vasco, Apdo. 644, 48080 Bilbao, Spain 
2 Departamento de Fisica Bdsica, Universidad de La Laguna, Spain 

The transient response of a stationary state of a quantum particle in a step potential to an 
instantaneous change in the step height (a simplified model for a sudden bias switch in an electronic 
semiconductor device) is solved exactly by means of a semianalytical expression. The characteristic 
times for the transient process up to the new stationary state are identified. A comparison is made 
between the exact results and an approximate method. 
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I. INTRODUCTION 

Assume that due to a constant flux of incident quan- 
tum particles a stationary scattering state is formed for 
a given one dimensional potential profile, and that the 
asymptotic potential level is changed suddenly. A phys- 
ical realization would be an abrupt change in the bias 
voltage of an electronic device |. The wave will then 
respond to the potential switch until a new stationary 
state is attained for any finite position x. Obtaining 
the characteristic time(s) of the transient is clearly of 
practical interest to determine the transport properties 
of small mesoscopic structures, but modelling the pro- 
cess by means of a grid discretization of space in a "finite 
box" is far from simple Q. The problem is that the 
boundary conditions at the box edges are not known a 
priori and involve simultaneous injection to and absorp- 
tion from the simulation (box) region. Some approxi- 
mate ways to deal with the transients have been proposed 
jl], ||, ||, H but, surprisingly, no exact solution has been 
obtained up to now. Our aim in this paper is to work out 
an explicit and exact solution of the transition between 
stationary monochromatic waves due to an abrupt poten- 
tial switch. While the calculation is performed for sim- 
plicity for a step potential that changes the step height 
suddenly, other potential profiles, e.g., containing square 
single or double barriers could be treated similarly. Our 
results are in fact applicable to the outer regions of an ar- 
bitrary cut-off potential with different asymptotic levels 
by inserting the appropriate reflection and transmission 
amplitudes. The basic trick to find the exact solution is 
to implement the action of the evolution operator of the 
new Hamiltonian on the initial state using an integral 
expression in the complex momentum plane obtained by 
Hammer, Weber and Zidell ||. 



where Q(x) is the step function and Vb > 0, a stationary 
state incident from the left has the form 



Mx) = h- 1/2 



e iq 0X /h + i^( ?0 ) e -i9ox/K ! x < o 

,ipox/ft ™ ^» n ~ 



\ x > 



where po = (g 2 . + 2m Vb) 1 / 2 and qo are positive, and the 
reflection and transmission coefficients for left incidence 
are given by 



go - Po 
Qo+Po'' 
2qo 



qo +Pa 

If the potential changes suddenly to 

v(x) = -K,'e(x) 



(3) 
(4) 

(5) 



at time t = the wave function will evolve in time. Find- 
ing ip{x,t) is equivalent to solve separately the time evo- 
lution of the initial functions 

ip 1 {x,t = 0) = h- 1 / 2 e t i° x / h e(-x), (6) 
V> 2 (M = 0) = h^^e^'^^&i-x), (7) 
i/)s(x,t = 0) = h^^e^^Qix), (8) 

and combine linearly the results with the appropriate co- 
efficients: ?p(x,t) — tpi(x, t) + R l ip2{x, t) + TQips(x,t). 

Of course the initial state could be different, in partic- 
ular a state incident from the right. Clearly, to treat any 
possible initial stationary state we have to consider also 
a fourth truncated plane wave: 



M^,t = 0) = h^^e-^^'Qix). 



(9) 



II. OBTENTION OF THE EXACT EXPRESSION 



For a potential step of the form 



V(x) = -V G(x), 



(1) 



Moreover, we should also allow for the possibility of a 
purely imaginary qo in ip2 to describe the evolution of 
an initially evanescent wave. We shall now examine the 
four possible cases. Note that the evolution of each of 
these initial states is a realization of Moshinsky's shutter 
problem [|| for the step potential. 



2 



A. ^i: initially a positive-momentum cutoff plane 
wave in x < 0. 



B. tp2' initially a negative- momentum cutoff plane 
wave in x < 0, or an evanescent wave. 



by 



The momentum representation of ipi(x, t = 0) is given The momentum representation of -02 (x, £ = 0) is given 



ii(g,t = 0) = 



i 1 
27r g — qo + iO 



(10) 



Following H , we shall rewrite the eigenstates of the new 
Hamiltonian, 

_ 1/2 f e ^+^(g)e-^/ n , x < r , 
^ _ h \ e**/*T l (q), x>0 ' [LL) 



in the form 



where 



^ q (x) =h- 1/2 e iqx/n +Ki, 



(12) 



by 



i(«,t = 0) = 



1 



2n q + qo + i0 



(20) 



with the pole again in the lower half g-plane. Following 
the same procedure used for 0>i and using the same 1Z- 
function decomposition [Eqs. ( |l2] ) and (fl3|)], the eigen- 
functions ip q , and the same contour contour, Ci = C\ 1 
one obtains 



C-2 



dqip q (x)4> 2 (q,t = 0) 



(21) 



In the evanescent case, E q < and qo = i(2mVo — p 2 ,) 1 / 2 , 
the pole lies in the lower imaginary axis. 



R _ h -i/2( R'e-^\ x<0 

X ' ~ 1 J<l £ ipx/H _ e iqx/h^ x > q l 1J J 



C. ^3: initially a positive momentum cutoff plane 
wave in x > 0. 



and 



E, 



= q 2 /2m = p 2 /2m-V \ 
= [p 2 -2 m Vtf' 2 , 
= [q 2 + 2mV f/ 2 . 



(14) 
(15) 
(16) 



The zero of energy is set by convention at the left level, 
and the amplitudes R l and T l are given in Appendix A. 
The square root in the definition of q is chosen with a 
branch cut in the p-plane between the branch points p = 
±(2?7i1q') 1 / 2 , whereas p has a branch cut in the g-plane 
between q = ±i{2mV ') 1/2 . This means in particular that 
q and p have the same sign for E q > and — V Q ' < E q . 
The solution method is based on writting the initial state 
as 



The treatment of "03 (x) is different. The momentum 
representation of 03 (x, t = 0) is given by 



b 3 (p,t = 0) = 



-i 1 



2ir p — po — iO ' 



(22) 



with a pole in the upper half-plane. We shall use the 
eigenstates 



/2 / T r (-p)e iqx / n , x < 

Jpx/h + Rr(_ p j e -i v x/H t x > 



(23) 



where the amplitudes are given in Appendix A. Similarly 
to Eq. (|l^) we write 



p (x) =h- l / 2 e ipx / n +1l 3 , 



(24) 



^l(x,t = 0)= / dqi/> q (x)M<l,t = 0), ( 17 ) 
JCi 

where C\ goes from —00 to +00 above all singularities 
(branch cut and pole). This is possible because 



with 



f dqMl,t = 0)TZx = 0, 



(18) 



as can be seen by closing the integration contour with a 
larg arc in the upper (/-plane and using Cauchy's theorem. 
Since q is an eigenstate of the Hamiltonian (even for 
complex q), the time dependent wavefunction is given by 

Mx,t)= f dqi> q (x)<f>(q,t = 0)6-^ ' h . (19) 



_ 1.-1/2 / T r (-p)e iqx / R - e ipx / n , x<0 



R '' " '' R r (-p)e-*P x / R , x>0 

The integral 



dpfo(p,t = 0)^3 = 



c 3 



(25) 



(26) 



vanishes for C3 going from —00 to 00 passing below the 
singularities (pole and branch cut). Note that the con- 
tour must now be closed in the lower half plane to apply 
Cauchy's theorem. Finally, 

^{x,t)= f dpMP,t = 0)M x )e~ lEqt/h - (27) 
Jc 3 
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D. ^4: initially a negative momentum cutoff plane 
wave in x > 0. 

The treatment is essentially the same as for ^3, but 
with 



where G(u) = g[k(u)]. This function has a simple pole at 
uq = (fco + b/2a)/f and possibly a branch cut, whereas 
ko is given in Table 1. It is now useful to separate the 
pole and branch-cut contributions explicitly and write G 

as 



Mp,i = 0) 



(28) 



G(u) = 



Ap/f 

U — Uq 



H(u), 



(34) 



III. CONTOUR DEFORMATIONS 



where Ao/f is the residue of G(u) at u — uq and the 
remainder, H(u), is obtained by substraction. Note that 
H{u) is either an entire function, if there is no branch 
cut, or its only singularity is the branch cut. 

The integral X is thus separated into two integrals, 
X = X' + X" . The first one may be reduced to a known 
function by deforming the contour along the steepest de- 
scent path (real-it axis) and taking proper care of the 
pole contribution, 



The explicit expressions for ijjj{x,t), j = 1,..,4 are 
given in Appendix B. The different terms contain expo- 
nentials of the form exp(±iqx/h) or exp(±zpa;/?i) (for 
terms with support x < or x > respectively). In each 
case the integral is better solved in the corresponding 
plane, q or p, by contour deformation along the steepest 
descent path to be described below. We shall generically 
use the variable k in both cases. Note that due to the de- I = 

composition of the stationary wave functions into three 
terms (associated with a T-amplitude, an i?-amplitude 
and an independent term, /, see Eqs. ([ll]) and (^3)) each 
wave function %pj(x,t) may be separated into three con- 
tributions that we shall denote as tpj a > where a — I,T, R: where w(z) = exp(— z 2 )erfc(- 



2 / Rt A 



du ■ 



Uq 



-me imx2 / M A w{-u Q ), 3 = 1,2 
iire imx2 / nt A w(u ), j = 3,4 



(35) 



ipj = ^2 a tpj.a- There are twelve of these terms, each with 
support in one half-line, and therefore twelve different in- 
tegrals. We shall denote as kj the poles in the momentum 
representation of the initial state ipj(x,t = 0). They are 
listed in Tabic 1 together with many other features of the 
twelve terms. The twelve terms may be written in the 
compact form 



-iz). In general the second 
integral, which involves the remainder H 7 has to be eval- 
uated numerically. 



1" = e l 



' m f / duH(u). 



(36) 



(29) 



dke-^ k2+kb ^g ja (k), (30) 

and vanish outside their support region. In the above 
expressions C1.2 = ij2i:h v l' 1 , 03,4 = —i/2ith}/ 2 , a — 
t/(2mh), 



However, the computational effort is greatly reduced by 
deforming the contour along the steepest descent path 
too. The branch cut, whenever it is present, cannot be 
crossed and has to be surrounded. This occurs for 



\mx/t\ < {2mVo) 1/2 . 



(37) 



Otherwise there is no branch cut contribution and the 
integral may be expressed as a series by expanding H(u) 
around the origin and integrating term by term, 



F = Q{x)(e ltv o /h -1) + 1, 



(31) 



I" 



Jmx /ht f^ 1 / 2 



and g(k) has always a pole (we shall drop the subscripts 
j,a unless they are strictly necessary). Moreover, the 
R and T-terms have also a branch cut singularity. The 
saddle point of the exponent is at k — —b/2a (xm/t for 
/ and T-terms and —xm/t for i?-terms) and the steepest 
descent path is the straight line Im(fc) = — (Re(fc)+b/2a). 
By completing the square, introducing the new variable 



H(u = 0) 



(38) 



+ 



OO 

E 

71=1 



1 x 3 x ... x (2n - 1) 
2 n {2n) 



H {2n \u = 0) 



u, 



u=(k + b/2a)/f, f = {l-i){mn/t) 1 ' 2 , 



(32) 



which is real on the steepest descent path and zero at the 
saddle point, and mapping the contour to the u-plane, 
the integral takes the form 



In practice the first term gives already a very good ap- 
proximation, even if Eq. ([57]) holds. The ^''-integrals 
have to be calculated numerically only for \mx/t\ < 
(2toV / / ) 1 / 2 and their relative importance with respect to 
w-terms from X' is only significant for rather small x at 
intermediate times, since as t — > 00, X" — > in all cases. 
A general analytical approximation making use of X' and 
the first term in X" is given by 



X = e} 



'/ht 



du t 



G(u), 



(33) 



Fe l 



2 /2ht 



f 

±Aqw(tu ) + ~^H{u = 0) 



(39) 
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FIG. 1: Exact density versus t (solid line), and approximation 
using Eq. ( p9| ) (dotted line) for x — lOOnm. Also shown is 
the contribution of X" multiplied by a factor of 10 (dashed 
line), and the point when the steepest descent path crosses 
the branch point (triangle). Mass=0.067m e , E q — 0.3eV, 
V = 0.3eV, V ' = 0.8eV. 



with the upper sign for j = 1,2 and the lower sign for 
j = 3, 4, see Fig. El Note the basic role of the w- 
functions, which may be considered elementary transient 
mode propagators of the Schrodinger equation || . They 
will show approximate wave fronts when x(uq = 0) lies 
within the domain of the term (i.e., when the saddle 
meets the pole). This occurs (for q n > 0) for the terms 
IT, 1R, 21, 31, 4T and 4R, see two examples in Figure [|. 

The long time behaviour may be obtained from the 
asymptotic (large-z) formula, 



w(z) 



2e~ 



Imz > 0, 
Im z < 0. 



(40) 



All -02.3 terms vanish as t — > oo for finite x, since these 
waves move initially away from the origin. In spite of 
this dominant motion, note that there is a transitory and 
generally small contribution of ip2T at positive x and of 
at negative x. On the contrary the w-functions of 
-01,4 pick up the exponential contribution in Eq. ( |40| ) 
which gives the new stationary states. In particular, as 
t — ► oo and for finite x, 



i^i(x,t) 

1p4,(x,t) 



e -iE q t/h 

e -iE q t/h 

h 1 / 2 



e iq x/H + R^ qo ) e ^iox/h^ X<Q 
T l (q )e l Po x / h 1 x>0 

T r (p' )e- iqoX / n , x<0 
e -i P ' x/n + Rr(p' Q yp^/r^ x > 



IV. EXAMPLES 



Figure || shows a typical wavefunction "density" ||2^| 
versus x at two fixed instants ti > t± > for Vq > Vq > 0. 
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FIG. 2: \ipi T \ 2 , \ip2Tf and |^i 3 /| at i = 10 fs versus x. 
mass=0.067m e , E q = 0.3eV, V = 0.3eV, Vq = 0.8eV. 
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FIG. 3: 



x (nm) 



versus x for t — 10 fs (upper figure) and 



t = 50 fs (lower figure). Exact solution (solid line). The 
diamond marks pot/m; the circle p' t/m; the square —qot/m. 
mass=0.067m e , E q = 0.3eV, V = 0.3eV, Vq = 0.8eV. 



For x > the main features are two flatter regions rep- 
resenting the old (to the right) and new (to the left) sta- 
tionary regimes separated by an oscillating structure. A 
simple semiclassical picture provides a good zeroth order 
explanation: assume a stationary flux of classical parti- 
cles in the old potential, with incident momentum qo and 
transmitted momentum pq. After the potential switch at 
t — 0, the last transmitted particle with momentum po 
will be at pot/m, whereas the first transmitted particle 
with momentum p' will be at p' t/m. Since p' > pq 
there is a region of width (p' — po)t/m where the two 
types of particles coexist. In the corresponding quan- 
tum scenario one may expect interference and oscilla- 
tions of wavelenght 2Trh/(p' — po) in this region, whereas 
the regions dominated by only one plane wave the den- 
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FIG. 4: h x Uav versus time at x — lOOnm, where the av- 
erage local frequency is defined as uj a v = —lm[(dtp/dt)/tjj]. 
Mass=0.067m e , E q = 0.3eV, Vo = 0.3eV, V ' = 0.8eV. The 
diamond and circle indicate xm/p' and xm/po respectively. 
The initial value is Vo + E q — Vq whereas the final asymptotic 
value is E„. 



sity does not oscillate and is proportional to the corre- 
sponding transmission probability, either |To(g )| 2 for the 
"old wave" or |T'(go)| 2 for the "new wave", which in the 
present case is smaller than the former, a somewhat sur- 
prising feature of quantum scattering off potential steps 
from the perspective of classical mechanics. The two cry- 
tical positions are marked in the figure with a diamond 
and a circle. The average local frequency Jl^, [n] shows 
also the transition between the two regimes, see Fig. ^. 

The x < region is clearly divided into "old" (to 
the left) and "new" regimes where incident and reflected 
components interfere. They admit also a simple anal- 
ysis: since the reflected wave stays dominated by mo- 
mentum — qo, the interference pattern wavclcnght stays 
the same in the new and old regimes, (nh/qo), and the 
only difference is the amplitude change due to the change 
of reflection probability from \R l (q )\ 2 to |i?'(q )| 2 - The 
transition at x = —q t/m is marked with a square in Fig. 

I 

In the example shown the dominant terms for x > 
are i/j^t and ipu representing respectively the "old" and 
"new" wave. Eq. ( |39| ) with these two terms only pro- 
vides a very good approximation. For x < the dom- 
inant terms are (incident wave), ip 2 i (old reflected 
wave) and ipm (new reflected wave). Again, the analyti- 
cal approximation describes the main features correctly. 
We may expect a worse performance of the analytical ap- 
proximation in processes involving tunnelling or evanes- 
cent waves, with the pole lying close to the branch cut: 
for example, when E q < 0, or —Vq > E q > 0. Some of 
these processes and the corresponding time scales have 
been studied recently in |l2| so we shall not insist on 
them here. 

Fig. |^ shows the density for a case in which Vq > Vq > 



FIG. 5: \ip\ 2 versus x: mass=0.067m e , E q = 0.3eV, Vo = 
0.8eV, Vq = 0.2eV. The square is at —qot/m, the circle at 
p' t/m and the diamond at pot/m. 
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FIG. 6: J x 100 at 
mass=0.042m e , E„ — 



the left box edge, x 
0.04eV, Vo = 0.42eV, 



= -22.48nm. 
Vo' = 0.62eV; 



box lenght: L — 44.96nm; number of grid points: N = 10 ; 
time step: At = 10 _6 fs. The circle is at t = xm/qo. The two 
stright lines indicate the values of J x 100 for the initial and 
final stationary functions. 



0. The new wave moves now at a slower pace than the 
old one so there is no interference structure between the 
two. 



V. COMPARISON WITH AN APPROXIMATE 
METHOD 

As stated in the Introduction, the basic difficulty to 
deal with transient phenomena between stationary scat- 
tering states by means of grid methods is that the time 
dependent boundary conditions at the box edges are not 
known a priori. Several approximate schemes have been 
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FIG. 7: J x 100 at x = 0: exact (solid line) and for two 
different boxlenghts: L = 179.84nm, iV=80000, At = 5.e-5fs 
(dashed line); and L — 44.96nm with the same grid density 
and At as in the other box (dots). mass=0.042m e , E q = 
0.04eV, V = 0.42eV, Vq = 0.62eV 
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FIG. 8: J x 100 at the right box edge, x = 22.48nm. 
mass=0.042m e , E q = 0.04eV, Vo = 0.42eV, Vq = 0.62eV; 
box lenght: L = 44.96nm; number of grid points: N = 10 4 ; 
time step: At = 10 _6 fs. The circle is at xm/p' and the dia- 
mond at xm/pQ. The two stright lines indicate the values for 
initial and final stationary functions. 



Substituting this form in the Schrodinger equation and 
neglecting the second order derivative of B, 



at 2m to ox 



(42) 



The first derivative of B is then calculated numerically 
at each time step with the first two spatial points, and is 
used to update the boundary condition for the next time 
step as 



iE qo At/n + go dB(x, t) c -ig x/n^ t ^ 

TO dx 

(43) 

Similarly, the wave function at the right edge grid points 
is written as 



Vv = C(x,t)e lp '° x/h . 



(44) 



Assuming again that C is linear in x and evaluating its 
derivative numerically with the two last grid points the 
boundary condition to the right is updated as 

+ At) r » yj(t) r e~* E ™ At/R - Po dC [ x ^ e -ip'o^At. 

m ox 

_ (45) 

We have adapted Koonin's grid method |l3| based on 
the Caley transform to this boundary-condition scheme 
and have calculated the "flux" |26| versus time at the box 
edges (Figures ^ and ||) and at the center x = (Figure |^) 
for the same potential jump considered in pi; the effective 
mass is taken as ttt, = 0.042 au (for Ino.53Gao.47As-AlAs) 
and the inicident energy corresponds to the Fermi level. 

The comparison with the exact results demonstrates 
that the linear ansatz is quite good at the right edge but 
fails at the left edge, where incident and reflected com- 
ponents interfere. The error introduced at the left edge 
propagates and affects eventually to whole simulation do- 
main, in particular the flux at the origin is deformed with 
respect to the exact one quite rapidly. Enlarging the box 
retardates the deviation from the exact result, see Fig. 
[j], but the computational cost becomes exceedingly large 
to reproduce correctly the whole transient at the origin. 



VI. DISCUSSION AND CONCLUSIONS 



proposed to overcome this difficulty, and our exact so- 
lution provides a needed reference to test their valid- 
ity and/or range of applicability. We have in particular 
made a comparison with a method proposed by Mains 
and Haddad based on looking a short distance into 
the simulation domain to determine what is coming out. 
Especifically, at the left edge region the wave is written 
as 

%ln = Ae lqox/h + B{x, t)e- iqoX ' h , (41) 



Since its discovery by Tsu and Esaki, tunneling 
through semiconductor nanostructures has been the ob- 
ject of a great attention due to its possible applications 
to ultrahigh speed electronic devices With the de- 
velopment of novel semiconductor nanostructures, it has 
become important to carry out theoretical and experi- 
mental studies on the tunneling process of carriers when 
an external bias is applied. In this way, electric field- 
induced electron transport has been recently explored in 
quantum dot arrays p2| , resonant tunneling diodes p3[ , 
and semiconductor superlattices ]24| . However, we note 
that one remaining key question in these experiments is 
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the analysis of the device transient response to an in- 
stantaneous potential step switching. The characteristic 
time of the response is of practical interest to determine 
the nanostructure transport properties and its possible 
applications to novel ultrahigh speed semiconductor de- 
vices. 

We have obtained an exact solution of the transition 
between two stationary scattering states due to the sud- 
den change in a potential step. Equivalently, we have 
solved exactly the Moshinski shutter problem for an arbi- 
trary cut-off plane wave in the step potential. (For other 
potential shapes see |l4|, [l6[ [l7, 0). The explicit ex- 
pressions obtained for their time evolution would in fact 
be directly applicably to an arbitrary cut-off potential 
with different asymptotic levels by using the appropriate 
transmission and reflection amplitudes. (For recent work 
on step-like potentials scattering see 19, 2(J). The exact 
results allow to identify characteristic times for the tran- 
sients. They also provide a needed reference for testing 
approximate methods that model time dependent open 
systems (finite systems exchanging particles with the out- 
side) with injecting and absorbing boundary conditions. 



for positive values of the arguments. The analytical con- 
tinuations for negative arguments are the amplitudes for 
the "outgoing" or "time-reversed" stationary states. 



APPENDIX B: WAVE FUNCTIONS 



These are the time dependent wave functions corre- 
sponding to the initial conditions given in Eqs. ([||^) for 
the potential of Eq. (ph : 
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iE q t/h /h q-p iqx /h\ < n 



2-K\fh 1 J c _dp 



(q+p)(qTqo+iO) 



x>0 



(Bl) 

where the minus and plus signs correspond to ip\ and ip2 
respectively, and 



APPENDIX A: REFLECTION AND 
TRANSMISSION AMPLITUDES 



The reflection and transmission amplitudes in station- 
ary waves with left and right incidence are given by 



R r {p) = T r ( P ) 2b- 



(Al) 



^3, 



2qe 



qx/h -iEgt/h 



(q+p)(pTPo-io) ' 

iEqt/h 



x<0 



2nVh \ f 0j d pg^(e**l* + f^e-W*), x > 

(B2) 

with the minus sign for ips and the plus sign for ip^. 



q+p 
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TABLE I: Features of the terms ipj. 



